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Abstract 



Modern experiment has shown that bulk matter exhibits quantum struc- 
ture and that this structure is electromagnetic in character. Therefore, we 
treat the electromagnetic field as the source of gravitation. In this way, the 
theory of gravitation becomes consistent with the quantum theory of mat- 
ter, which holds that electric charge (or 'generalized charge') is the most 
fundamental attribute of matter. 

The following predictions of the theory may lie within reach of present- 
day methods in astronomy: 

(1) any massive body generates a time-dependent gravitational field; 

(2) there is a linear correlation between the gravitational red-shift of a 
stellar source and the energy of cosmic rays emitted by that source, 
given by Av/i^q = energy (eV)/10 27 ; 

(3) the maximum energy of cosmic rays is 10 27 eV; 

(4) this limit is associated with an infinitely red-shifted stellar object — an 
"electrostatic black-hole" — at the potential c 2 /G 1 / 2 = 10 27 volts. 

Finally, the theory predicts that the gravitational potential near any charged 
elementary particle is many orders of magnitude greater than the Newtonian 
value. 
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I. Introduction 



In the theory which follows, electromagnetism is taken to be the source 
of gravitation. The foundation for this hypothesis lies in the observed na- 
ture of matter itself. Modern experiment has shown that matter can no 
longer be described, in any fundamental way, by an amorphous 'density of 
mass.' Rather, the experiments reveal an intricate quantum structure which 
is largely electromagnetic in character. This profound change in our knowl- 
edge of matter has met with no corresponding change in the theory of the 
gravitational field. Consequently, a perplexing estrangement has replaced 
the close historical tie between gravity and matter. It is toward restor- 
ing this close tie that we propose to trace the origin of gravitation to the 
quantum of electric charge. 

A development of this nature cannot take place within the context of 
Einstein's theory of general relativity. According to that theory, electro- 
static energy produces a repulsive gravitational field. (This follows from the 
Reissner- Nordstrom solution of Einstein's field equations, for a free charged 
particle, at rest.) In order to arrive at the attractive gravitational field of 
experience, we must reformulate the problem of space and time. To this end, 
we abandon the four-dimensional vector basis of general relativity. At each 
point in space-time, we introduce a scalar basis eo and a three-dimensional 
vector basis e$, i = 1,2,3. We begin by expressing the Lorentz transforma- 
tion in terms of this scalar, three- vector basis. This was first accomplished by 
L. Silberstein in 1912 [1]. We then go on to eliminate relative motion from 
the theory, by suitably restricting the transformations of space and time. 
This gives rise to a covariant expression which is identified with the grav- 
itational field. Thus, according to our theory, a gravitational field induces 
non-uniformity in space and time which 'breaks the symmetry' of Lorentz 
invariance. This field has well-defined energy, momentum, and stress at each 
point in space. Moreover, it has tensor character and cannot be arbitrarily 
transformed to zero. Finally, it interacts with electromagnetism via a tensor 
force and thereby justifies our unified treatment of the field. 



1 Evidence for this may be found in any book on modern physics. On the one 
hand, there are the strong, electromagnetic, and weak interactions between the 
elementary particles of matter; on the other hand, some 40 orders of magnitude 
removed, there is the Newtonian gravitational interaction of 'bulk' matter. 
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II. Transformations of Space and Time 

The quaternion formulation of special relativity was introduced by L. 
Silberstein in 1912-1913 [1]. In this work, Silberstein made use of the original 
quaternion basis due to Hamilton (l,i,j,k), where ij = — ji = k (cyclic) 
and i 2 = j 2 = k 2 = —1. Here, we adopt the system of Pauli operators 
(cjo, (T±, <T2, <r 3 ), where <ti<X2 = — <x 2f \ = icr 3 (cyclic) and a\ = cr\ = n\ = 
o"q- The fundamental interval in this theory is a quaternion, i.e., a sum of 
scalar and three-vector displacements 



ds = cdr + dr 

= (Jo dx° + cr i dx' 1 

= <T„dx» (1) 

The operator uo commutes with all operators and, in this way, time is 
distinguished as a scalar variable. The interval ds transforms according to 
the general rule 

ds = a(ds)b (2) 
where a and b are also quaternions. For example, the Lorentz transformation 

= j^-^dx 1 ) (3) 

dxt = 1 (dx 1 --dx°) (4) 
c 

dx 2 = dx 2 (5) 
dx^ = dx 3 (6) 

is given by (2), with [1,4] 



a = b 



= ±-i=(cr ov /7 + 1-0-1^/7-1) 



±(o- cosh| 



(j\ sinh 

T 



(7) 



and 
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-1/2 



7 = cosh£ = 1 - — 



A scalar product can be defined by introducing the conjugate basis 
(<7o, — <7j) and conjugate interval [5] 

ds = dx^ = oq dx° — <Ji dx i 
It follows from the commutation relations that 



dsds = (oodx + o-idx l )(oodx — o-j dx^) 
= oo(dx° — dx 1 — dx 2 — dx 3 ) 

Transformers a and b have unit magnitude 

aa = bb = o"o 
therefore, the scalar product is invariant 



ds ds = a(ds)b b(ds)a = ds ds 
It is convenient to form the symmetric products 

1 



where 



Then 



Vnu 



( 1 



\ 



ds ds = (JpOy dx^dx u 

= 2 + °V0/J dx^dx u 

= o"o tyx" dx^dx v 



The above formalism may be extended in a straightforward manner to 
variable basis systems: 



ds = cdr + dr 

= cq(x) dx° + Gi{x) dx 1 
= e M (x) dx^ 



(16) 



Here, we retain the scalar nature of time and the three-vector nature of 
space by setting 



eo(x) = e° (x)<Jo 



(17) 



ei(x) = e j i(x)o-j 



The conjugate interval is 



ds = e M dx^ = eo dx° — dx 1 
Symmetric products are found to yield scalar functions 



where 



and 



1, _ 



f goo o o o \ 



gij 

V o J 



goo 

9ij 



e°oe°o 



These metric coefficients yield the invariant product 



(18) 
(19) 
(20) 



(21) 



(22) 
(23) 



ds ds = e^ev dx^dx v 
= cr g^u dx^dx v 



(24) 



2 The term "metric" will be used, even though component goo has no foundation 
in geometry. 
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III. The Structure of Space and Time 

The basis system e^(x) varies from point to point in a given manifold, 
the vectors ej changing in both magnitude and direction. This behavior is 
expressed in terms of coefficients [] Q^ x 

= exQi* (25) 
The rate of change of the scalar basis must be of scalar form 

Vue = e Ql (i/ = 0,l,2,3) (26) 

Qi„ = (i = l,2,3) (27) 
while the rate of change of the vector basis is a three- vector 

V^e, = ej Ql (28) 

Qi = (29) 

Therefore, the are generally non-symmetric in lower indices. 

Certain of the coefficients Q^ x have tensor character. These are deter- 
mined in the following way. A given observer may record a time measure- 
ment dr and vector displacement dr with respect to basis e M : 

cdr = eo dx° , dr = e$ dx l (30) 
This observer is free to change the vector basis — ► e/, in which case 

dr' = e/ dx 1 ' where dr' = dr (31) 

He may also change the scalar basis eo — * eo' (rate of time coordinate 
clocks), without affecting the length of time involved 

c dr' = eo' dx ' where dr' = dr (32) 

The new time coordinate is independent of space coordinate labels, x ' = 
x° (x°), and the new coordinates are independent of clock rates, x l = 

3 In the literature, Cartan's notation is often used, de M = e\Q^ u dx v , even though 
the e M are not generally functions and the de^ are not differentials. Here, we set 
de^ = (V v e ll )dx v so that (25) follows. 
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x i (cc- 7 ). In short, since no relative motion is involved, the transformation 
reduces to 

dx° 

e °' = dx^ e ° (33) 

dx^ 

e/ = ^ (34) 
Applying these restricted transformations to the we find 

/ dx^ dx^^ x^ 

QOO = d^ Q ^ + -^0^2 ( 35 ) 

Q°ok = ^Ql (36) 
ox K 

■ i dx^ dx^ dx^ 

Q 5° = dx^dx^dx^^ (37) 



dx 1 ' dx b dx c „„ dx 1 ' d 2 x a 



Ql]k dx- dxo' dx k ' Qtc + dx- dxi'dx*' (38) 

Functions Q^ k and Q* transform as tensors and thus have invariant signif- 
icance for our observer. Defining 

we obtain the complete array of tensor components 

(40) 
(41) 
(42) 
(43) 

The vectors e, form the basis of a time-dependent, three-dimensional 
geometry, which is characterized by metric components gij. The metric 
changes from point to point in the spatial manifold according to the formula 

^ = 9iaQ% + g ja Qt k (44) 

By imposing nine symmetry conditions 



Q[0k] 


— Qok 


— Qlo 


— Qok 


Q°[jk] 


= Q% 


— Qkj 


= 




= Q)o 


— Qhj 


= Q}o 


Qyk] 




- Qlj 
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Q)k = Qlj (45) 
this formula can be inverted in the usual way to yield 

The time-dependence of the metric is expressed in terms of coefficients Q'j : 

§|=ffiaQ?0+*aQS) (47) 

In order to invert this formula, we impose three conditions 

9iaQ° = QjaQl (48) 



obtaining 



<& = 5«-^f («) 



Finally, the function qqq changes according to 



~ ^Ql, (50) 



which immediately gives 



O o _ I oo^oo r _ n 

^" 2 5 ^ (51) 

These formulae express all 28 independent coefficients in terms of the 
28 derivatives dg^/dx^. 

We can state these results more concisely by noting that conditions (48) , 
together with (45), are equivalent to the equation 

a^QfvX] + 9vaQ[\n\ + 9XaQ^ v] = o (52) 
We now begin with the general formula 

-^x = 9^aQZ\ + 9uaQ^\ (53) 

in the combination 
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V dx x dx u dx^ J 

= 9v P Q^x + g^Qtx + g^pQxu + 9\ P Q% - 9\ P Q p V p - 9u P Q p Xll 
= 2 9 P pQ P vX + 9p P Q P [Xv] + 9v P Q p m + 9\ P Q P [M 

and then make use of (52) to find 



ivp, 



+ 



dx x dx u dx^ J 
Raise index /j, and define the Christoffel symbols 

= 1 na ( d 9ua , dg a x 

vX 2 9 \dx x dx v dx a 

in order to obtain the formula 

Kx = Qlx + 9^9xpQ[ 



[H 



au] 



IV. Conservation of Energy-Momentum 

The density of energy, momentum, and stress is represented by 



T = e M <g> e u T^ 



(54) 



(55) 



(56) 



(57) 



(58) 



where T» v = T Ufl . Consider the product 



(V=jT, dV) = e^gT^dV u 

= e V=g~ T 00 dV + e V=9 T 0j dV 3 
+ ej v 73 ^ T i0 dV + ei V^g T ij dVj 



(59) 



The first term in this expression is a scalar, which gives the amount of energy 
in region dVo 



dE = e ^g~T m dV 
while the momentum is given by the three- vector 

cdP = ei^g~T i0 dV 



(60) 
(61) 
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The remaining terms represent the flow of energy and momentum through 
surfaces dVj. 

The divergence theorem for energy-momentum is [6] 



= / e 



1 dyf^gTV* 
/:r 3 dx v 



+ 



Energy-momentum is conserved if 



divT"" = W J_. + Q^ x T vX = 



ax- " u (63) 

The divergence takes on a more useful form, if we replace by means of 
(57): 



divT^ 



1 dy/^gTV" 



+ 9 Winamp 



(64) 



(A semicolon indicates the covariant derivative defined in terms of the Chris- 
toffel symbols.) 

The density of electromagnetic energy, momentum, and stress is given 

by 



In order to investigate conservation, we calculate 

U1V 1 e-m — 1 ;u e-m ^ U ^[ijape-m 

which, by virtue of Maxwell's equation, becomes 

alvl e-m — r a-J ^9 W [ V a] P e-m 



(65) 



(66) 



(67) 



Therefore, electromagnetic energy-momentum is not conserved in regions 
free of charged matter, due to the interaction term in • In this way, the 
gravitational field becomes physically manifest. 



Performing a similar calculation with the matter tensor 
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(68) 



we find the covariant derivative 



1 ;uM - PC U ^ v -tC U 



1 dy 



dx u 



(69) 



The second term is zero, if rest mass is conserved. Taken together with (67), 
we obtain the divergence 



div(T^ + Trj = P c?(u^ + T» x u»u^+F\r 

+ ^QU^M + ^e-J (70) 

The first expression, set equal to zero, gives the equation of motion for 
charged matter 

PC 2 {f^ + lt A + F^r = (71) 

and we are left with the gravitational interaction 

div (TZ + 3Tm) = STO^.] ( T PM + 3?e-m) (72) 

This interaction immediately unites electromagnetism and matter with 
gravitation. Moreover, the interaction has tensor character, which implies 
that the gravitational field must have real dynamical content (energy, mo- 
mentum, stress) in each region of space. This is established in the following 
section. 

V. Gravitational Field Equations 

In forming a Lagrangian suitable for gravitation, we begin with the in- 
variant expression 

(G = 6.67 x 10~ 8 cm 3 /g-sec 2 ). The coefficient is chosen such as to yield the 
gravitational energy tensor 
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Tgrav = ^ {^V^f, A ~ ^9^^ } (74) 
The non-zero components of Q^a] were determined in section III: 

Qh = <& = ± 9 M t? (75) 

Q'w = ^ = (76) 

Consider the special case of a weak gravitational field, 



5oo = 1 + 4^ ( 77 ) 

9ij = ( 78 ) 

where ip is the Newtonian potential. Substituting into (74), we obtain the 
weak- field limit of T^ av : 



TfL = ^G^) 2 (79) 
^ = i^{^^ m ( V ^)( V ^)-^ lJ (VV) 2 } (80) 

T^? at , is the positive definite energy density of the gravitational field. T^ rav 
is the stress tensor of Newtonian gravitation [7]. The stress is compressive 
along the gravitational lines of force, with tension acting between the lines 
of force. 

Returning to the general case, we introduce arbitrary variations of the 
seven fields g^ v = (g 00 , g^) and calculate 



A x 



6 1 8ttG { 



' C 1 „00 Ja „mb ^9am ^9bl , 1 J m 00 „00 ^900 $500 1 /— ,4 



4 y y y dx° dx° 4 y y y dx l dx m 1 
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= /{: 



1 d 



d 
dx° 

We have made use of these formulae 



+ 



fi—— 



(V=9 9 00 Q)o) + 9uSg lj V=g~d 4 x (81) 



$9nv = -g m gup $g al3 



(82) 
(83) 
(84) 



V dx x J dx x 

as well as integration by parts; variations Sg^ u have been set equal to zero at 
the limits of integration. The first term in (81) will yield Poisson's equation 



V 2 V> = 4vrGp 



(85) 



if we introduce the matter tensor (68). This is accomplished by means 
of the variation 

5 J - {pc 2 (g^un 1/2 - pc 2 } ^-gd i x= l -j M 8gT V=2 d 4 x 

= \j {T$ M gooSg 00 + T) M gM 3 } d 4 x (86) 

The final expression to be introduced is the energy tensor of electromag- 
netism [8] 

5 J - l - Fflv F^^d A x = \J T^ e _ m Sg^^gd 4 x 
= \f {TSe- m gooSg°° + Ti e _ m gu5g l i } ^~g d 4 x 



(87) 



where T^ m is given by (65). 

Gathering terms together and setting coefficients of 5g m and 5g %: > equal 
to zero, we arrive at the field equations of gravitation 



1 d_ 

AirG y/^ ~g dx l 



{^gg lm Ql m ) + T^ = Q 



(88) 
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d 



rpiiv re p resen t s the total energy, momentum, and stress 

1 1 grav ' M ' e—m 

The corresponding invariant equation is 







1 



8 



4ttG y/^g dx l 



(V=gQl + T = o 



where 



(89) 



(90) 



(91) 



VI. The Gravitational Field of an Electron 



(92) 



According to our hypothesis, it is the electromagnetic nature of matter 
which gives rise to gravitation. Thus, our central problem is to determine 
the field of a discrete charged particle, at rest. The field equations admit an 
exact solution for this case, of the form 







\ 



(93) 



\ -r 2 sin 2 6> J 

where v and A are functions of r which vanish at infinity. The electromag- 
netic energy tensor is [9] 



T, 



\iv e—m 



8vrr 4 



If? 



V 



o 



— e 







9 • 2 

r sin i 



\ 



(94) 



where q is the electric charge in ordinary c.g.s. units (stat-coulomb). The 
components of the gravitational energy tensor are found from (74) and (93) 
to be 
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\ 



(A prime indicates differentiation with respect to r). Finally, we have 



(95) 



y -r 2 sin 2 6> j 



d_ 

dx l 



= _I r 2 sin0e (-A)/2 (n + 2 v , + l 
Substitution into (88) and (89) yields two distinct equations 



. 2 .j . V 2 - Va' 

2 



c 4 r 4 



1 ./2 , Gg 2 A 



4 c 4 r 4 



e A = 



(96) 

(97) 
(98) 



Adding the two equations gives 



or 



v " + l v < _ \ v 'x' = 
r 2 



2 1 

d(lni/) + -dr = - dX 
r 2 



(99) 



Upon integration, we find an equation of the form (98). In order to arrive 
at a unique solution, we set 

A = -v (100) 

That is, we invoke a reciprocal relationship between time dilatation and 
length contraction along the direction of the field. Equation (98) now gives 



e v l 2 dv = ^ ^ dr 
err 



This equation integrates to 



G 1 / 2 c 
c 2 r 



(101) 



(102) 
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the constant of integration being chosen such that e v = e A = 1 at infinity. 
The components of the metric tensor are 



/ 



1 - 



G l / 2 t 
c 2 r 



V 



o 



— r 2 sin 2 i 



(103) 



/ 



The gravitational field of an electron is given by setting q = e = 4.8 x 
10~ 10 stat-coulomb in (103). Noting that 



G l/2 e/c 2 = 1.4 x 10- 34 cm 
the field may be represented by a gravitational potential 



(104) 



(105) 



This potential is some 21 orders of magnitude greater than the Newtonian 
value, —Gm/r. It acts in those regions of space which are subject to the 
electrostatic field of the electron. 

VII. An Electrostatic Red-Shift in Stars 



In the previous section, it was shown that a spherically symmetric elec- 
trostatic field produces a gravitational potential given by 



i, = -G 1 ' 2 - 



(106) 



where q is the absolute value of electric charge. This potential, in turn, will 
produce a red-shift of light given by 



Av _ _ G 1 / 2 q 

u c 2 



c 2 r 



(107) 



where v$ is the light frequency in the absence of gravitation, i.e., at r — ► oo. 
In this section, we discuss a plausible observation of this red-shift and its 
correlation with ordinary electrodynamic processes. 
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For purposes of our discussion, consider the following hypothetical stel- 
lar object, in which a spherically symmetric charge separation has occurred. 
Neutrons and protons, held strongly in the core, are surrounded by elec- 
trons which have been forced away by thermal energy and mutual repulsion. 
Under the simplest assumption, all electrons lie in a spherical shell of fixed 
radius. This model may be treated in terms of our theory as follows. Each 
neutron produces a Newtonian potential 

Gm 10~ 31 , , 

Y>n = = 108 

r r 

(m = 1.7 x 10 _24 g), while each proton, according to (106) produces the 
much greater potential 

^ G i/2^ = _1P_^! (log) 

r r 

Ignoring the neutron contribution, we write the gravitational potential (106) 
in the form 

= -G 1/2 <f> (110) 

where <p is the total electrostatic potential. A ray of light, being emitted 
from the core and passing through potential (f>, will suffer a red-shift (107) 
given by 



= — ~ ( 

u c A 



(111) 



This linear relation holds true over a very wide range of electrostatic poten- 
tial (see below). In addition to light, one would expect positively charged 
particles (cosmic rays) to be emitted from such a stellar object. Particle 
energy is commensurate with the electrostatic potential. Thus, according to 
(111), one should observe a linear correlation between the red-shift of light 
and the energy of cosmic rays. 

The exact solution to the electrostatic problem (103) takes the form 
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9tiu 



G l/2 



V o 

Component goo gives the frequency relation 



-r 2 sin 2 6> J 



(112) 



— =%o 

"0 



-1/2 



G 1 ' 2 \ 1 



(113) 



and red-shift formula 



GV2 



Au _ c 2 



v » 1 G1/ \ 



(114) 



According to this formula, an infinite red-shift corresponds to a potential of 



c 2 /G 1/2 = 3.5 x 10 24 stat- volts 

= 10 27 volts (115) 

It follows that cosmic ray energy can be no greater than 10 27 eV. 

It would be of interest to search for evidence of a correlation between 
the gravitational red-shift of a discrete source and the energy of cosmic rays 
emitted by that source. This may be within reach of modern experimental 
methods in astronomy. Perhaps of greater significance is the predicted upper 
limit to cosmic ray energy, 10 27 eV. This limit is associated with an infinitely 
red-shifted stellar object, an "electrostatic black-hole," at a potential of 10 27 
volts. The astronomers have found the maximum energy of cosmic rays to 
be at least 10 20 eV. This compares favorably with the theoretical value. 

VIII. Concluding Remarks 

According to equation (89), any flow of momentum (i.e., stress) will give 
rise to a time-dependent gravitational field. This applies, in particular, to 
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the field generated by a stationary uncharged mass such as the Sun. In a first 
order approximation to this problem, all stress terms may be ignored; (88) 
then yields Poisson's equation and the static Newtonian potential. How- 
ever, in second order, the theory predicts the existence of time-dependent 
perturbations in the Newtonian field, due to the presence of material and 
gravitational stress. Recently discovered resonant oscillations in the Sun 
show that the solar field is, indeed, time-dependent. It remains to be seen 
whether this provides corroboration for our theory. 

The above behavior contrasts with that of the field surrounding a sta- 
tionary electric charge (section VI). Here, the electrostatic stresses are equal 
and opposite to those of the gravitational field. Therefore, the net flow of 
momentum is everywhere zero, and the solution is time- independent. This 
result underscores the significance of electric charge in our theory. Moreover, 
it calls to mind the central role played by electric charge and 'generalized 
charge' in elementary particle theory. The enormous strength of the elec- 
tron's gravitational potential may point toward an eventual link with the 
physics of elementary particles. 
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